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Abstract 

The AB s-valent dimer is used to analyse bond formation and charge transfer 
within the tight-binding (TB) approximation. In this way a physical interpretation of 
the electronic structure and binding energy within density functional theory (DFT) 
is obtained which lends itself to the derivation of covalent and ionic interatomic 
potentials. 
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1 Introduction 

The AB s-valent dimer provides the simplest system to illustrate the underlying ap- 
proximations and concepts behind the TB method. Whereas density functional the- 
ory (DFT) centres on the observable quantity p(r), the electron density at all points 
in space r, the TB method coarse-grains the problem in terms of chemically-intuitive 
but non-unique quantities such as the bond integral (3{R), and overlap integral 
S(R), defined in terms of the internuclear distance R between a given pair of atoms. 
In the following sections we will see how the TB method allows the DFT electronic 
structure and binding energy curves to be interpreted in terms of physical concepts 
such as covalency and ionicity. This leads naturally to a further coarse-graining of 
the problem in terms of covalent and ionic interatomic potentials. 



2 Electronic structure 

2.1 Energy levels 

Let us consider bringing together two s-valent atoms A and B to form the AB di- 
atomic molecule, as illustrated in Fig. 1 [1], This could represent, for example, either 
the homovalent dimer H2 (with Is valence orbitals interacting) or the heterovalent 
dimer LiH (with the 2s Li valence orbital interacting with the Is H orbital). The 
resultant electronic structure of the dimer is obtained by solving the corresponding 
Schrodinger equation 

- ^VVabM + V AB (r)iJAB(r) = E^j AB (r) , (1) 




Free Atom A Free Atom B AB diatomic molecule 



Figure 1: The potential of the AB diatomic molecule may be approximated by the over- 
lapping of the free atom A and B potentials. E° A and E B are the free-atom energy levels 
of the valence s electrons. The internuclear distance R runs from A to B. 



where we have used atomic units 

H = 1, m = 1, e = 1, 4tteq = 1 . (2) 

Thus, h 2 /2m = \ in Eq. (1) and the unit of energy is the Hartree (= 2 Rydberg = 
27.2eV) and the unit of length is the au = 0.529A. 

The TB method solves this differential equation by assuming a minimal basis of 
atomic orbitals, namely 

i>AB{r) = c A <p A {r) + c B (f>B(r - R) (3) 

where R is the internuclear separation with the coordinate system centred on atom 
A as in Fig. 1. The atomic orbitals satisfy the free-atom Schrodinger equations 

-\v 2 <t)A + V A <i)A = E A( t )A (4) 

-\v 2 <I>b + V B 4> B = E° B (i> B (5) 

with E\ and E B being the free-atom energy levels of the valence s electrons. The 
coefficients c A and cb can be found by pre-multiplying the Schrodinger equation for 
the dimer Eq. (1) by ipAB(r) , integrating over all space, and writing the resultant 
expression in the form 



/ 



i> AB (H - E) i, AB dr = , (6) 
where H is the Hamiltonian operator for the AB dimer, namely 

H = -^V 2 + V AB . (7) 
Then substituting Eq. (3) into Eq. (6) we recover the TB secular equation 

H AA -E H A b-ES \ ( c A 



Hba — ES Hbb — E J y cb ' 
where the on-site Hamiltonian matrix element H AA is given by 

H AA = J (f> A H(t> A dr , (9) 
and similarly for Hbb- The intersite Hamiltonian matrix is given by 

H AB = J (t> A Hct>Bdr, (10) 

and the overlap integral by 

S = J M B dr . (11) 
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The orthogonal TB model assumes that S = 0. Within Molecular Orbital theory Eq. 
(8) is referred to as the Hiickel equation for S = and the extended Hiickel equation 
for S / [2]. 

We will solve this TB secular equation for the eigenvalues and eigenvectors by work- 
ing with respect to the average on-site energy of the dimer, namely 



E=-{H BB + H AA ). 



Then, defining 



and 



Eq. (8) takes the form 



E-E 



Hbb — Haa , 



— e 



6-eS 

This has non-trivial solutions for 



(3-eS 



CA 
CB 



— e 



(3-eS 

where the bond integral j3 is given by 



(3-eS 



0, 



f3 = H AB -ES. 



(12) 

(13) 
(14) 

(15) 

(16) 
(17) 



We see that the absolute energy E with respect to the vacuum level does not enter 
Eq. (16) explicitly only the relative energy e = (E — E). The absolute energy E is 
hidden within the definition of /3; the physical importance of this grouping will be 
discussed in section 2.3. 



The exact eigenvalues of the TB determinantal Eq. (16) are given by 



e ± = 



-PS T V/3 2 + (l-S 2 )(A/2) 2 l /(l - S 2 ) . 



(18) 



In order to keep our analytic analysis of this TB model for the dimer tractable for 
A / 0, we now neglect all second and higher order terms in the overlap S since 
they only lead to third-order and higher contributions such as f3S 2 and AS 2 in the 
eigenspectra. The eigenvalues then take the form 



E ± = E + \p\S T^ 2 I (A/2) 2 , (19) 

since the bond integral j3 will be seen in section 2.3 to be negative. Thus, as illus- 
trated in the left-hand panel of Fig. 2, the eigenvalues comprise a bonding state 
and an anti-bonding state that are separated by 2|/3| and 2|/3ab| for the homovalent 
and heterovalent cases respectively, where (3ab = \/ ft 2 + (A/2) 2 - The shift in the 
energies due to the overlap repulsion \ fi\S is not shown explicitly in Fig. 2. 
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Figure 2: The bonding and antibonding states for (a) the homonuclear and (b) the het- 
eronuclear diatomic molecule. The shift in the energy levels due to overlap repulsion 
has not been shown. 

2.2 Charge density 

As illustrated schematically in Fig. 2, H2 and LiH have their two valence electrons 
occupying the bonding state iPab( t ) w ith anti-parallel spins. Substituting the eigen- 
values from Eq. (19) into the TB secular equation we find that the eigenfunctions 
are given by 



where 



V2 



l±(A-5)/Vl + A 2 



(20) 



(21) 



and 



B 



V2 



it(a + s)/Vi + a 2 



(22) 
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neglecting second and higher order contributions in the overlap. The normalized 
atomic energy-level mismatch A is defined by 

A = A/(2|/3|). (23) 

Thus, for a given internuclear separation R the number density of the valence elec- 
trons in the AB dimer can be written as 



PAB (r) = 2[^ AB (r)}\ (24) 

with the corresponding electronic-charge density = —ep AB (r) = ~ Pab (?) in atomic 
units. It follows from Eqs. (20) and (24) that 



PAB (r) = 2{c + A fp A {r) + 2{c + B fp B {r - R) + 4c+c+^(r)^ B (r - R) (25) 
= N APA (r) + N BPB (r - R) + 2Q AB (j) A {r)(j) B {r - R) , 

where 

PA(B)(r) = \4>A{B){r)] 2 ■ (26) 

N A and N B are the number of valence electrons on atoms A and B respectively, 
namely 

N A = 1 + A//l + A 2 - S/\l 1 + A 2 , (27) 

and 

N B = 1- A/7 1 + A 2 - 5/^ 1 + A 2 , (28) 

to first order in S. The prefactor Q AB in the last term of Eq. (25) is the bond order 
between atoms A and B, namely 

AB = l/^Jl + &- S. (29) 

Eq. (25) for the number density may be simplified by grouping the overlap contri- 
bution in the first two terms with the bond-order contribution in the last term to 
give 

PA B {r) = (1 + q)pA(r) + (1 - q) PB {r - R) + &p C ov(r) , (30) 

where 

Pcov {r) = 20 A (r)0 B (r - R) - S[p A (r) + p B (r)} . (31) 
It follows from Eqs. (27) and (28) that the charge q is given by 

q = A/^l + A 2 >0, (32) 
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which is positive due to our choice of Eb > Ea in Fig. 1. On the other hand, it 
follows from Eq. (29) that the bond order is given by 

9= 1/^1 + A 2 = ^l-q 2 . (33) 

We see from Eq. (32) that the normalized atomic energy-level mismatch A can be 
written in terms of the charge q as 

A = g/7TV. (34) 

Thus, A — > corresponds to the covalent limit with q = and 9 = 1, whereas 
A — > oo corresponds to the ionic limit with q = 1 and = 0. 

The nature of the covalent-bond density /? col ,(r) is illustrated in Fig. 3 for the hydro- 
gen molecule where we see, as expected, that the electrons flow from outside the 
bond region to inside. Importantly, the total covalent-bond charge is zero because 

J Pcov(r)dr = 2S - 2S = . (35) 

Interestingly, the net atomic charge of ±q is consistent with the definition of Mulliken 
charge where the overlap contribution 2@abS from Eq. (25) is assumed to be di- 
vided evenly between both atoms A and B [2]. It follows from Eqs. (27), (28) and 
(29) with the values of q and 6 from Eqs. (32) and (33) that to first order in S 

qMuiuken = (<? - ®S) + OS = q . (36) 

In this section we have demonstrated that a given choice of A and f3 in the TB secular 
equation leads to a specific value for the atomic charge q through Eq. (32). In the next 
section we will show that these TB parameters A and (3 are themselves functions of 
the charge so that Eq. (32) must be solved self-consistently for the charge q. 

2.3 Expressions for TB parameters 

The key parameters that enter the TB expression for the energy levels, namely the 
atomic energy-level mismatch A and the bond-integral ft, depend on the Hamilto- 
nian matrix elements Haa, Hbb arid Hab through Eqs. (14) and (17). These matrix 
elements in their turn depend on the potential Vab{t) through Eqs. (7), (9) and (10). 
Within DFT the potential seen by the electron in the Kohn-Sham equations can be 
written as the sum of three contributions 

V(r) = V ion (r) + V H {r) + V xc (r) , (37) 

where the first is the potential due to the ion cores, the second the Coulomb or 
Hartree potential from the valence charge density, and the third the exchange-correlation 
potential that enters within DFT. In this analytic treatment of the dimer we will ne- 
glect this latter contribution as it varies non-linearly with the density unlike the 
Coulomb term. We will assume that it can be subsumed into the TB parameters 
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Figure 3: The electron density of the homonuclear molecule (upper panel) can be re- 
garded as the sum of the non-interacting or frozen free-atom electron densities (lower 
panel) and the quantum mechanically induced covalent-bond density (middle panel). 
The dashed curve represents the first-order approximation, Eq. (31), for the bond den- 
sity, the deviation from the exact result (full curve) being due to the sizeable value of 
the overlap integral namely S = 0.59 at R = 2au. 

when they are fitted to the DFT eigenspectra and binding energy curves. We will 
see in section 3.3 that this is indeed a good approximation as it leads to a consis- 
tent physical picture for both the non-magnetic DFT energy levels and the binding- 
energy curves of the hydrogen molecule. 

Thus, using Eq. (30) for the valence density and neglecting the exchange-correlation 
contribution, the potential Vab{t) for the AB dimer takes the form 

V AB (r) = V A (r) + V B {r) - q[Vg(r) - Vftr)] + . (38) 

The first two terms are the Coulomb potential that results from overlapping neutral 
free-atomic potentials, as sketched on the right-hand side of Fig. 1. That is, 

V A (B)(r) = V^ B) (r) + V^ B) (r), (39) 

where Vi on is the potential for the ion core that falls off inversely with distance out- 
side the core region such that 
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V A T B) {r) = Z A{B) /\r-R A{B) \iox\r-R A{B) \ > R%(b) ■ ( 40 ) 

For our monovalent AB dimer Z A = Z B = 1 and Ra = 0,Rb = R from Fig. 
1. Whereas for hydrogen the core radius is zero, for lithium it takes some finite 
value enclosing the Is shell of core electrons. V p is the potential due to the valence 
electrons of the free atoms, namely 

V A(B)( r ) = J PA(B)(r' ~ R M B))/\r - Adr' . (41) 

The superscripts p in Eqs. (39) and (41) are used to remind us that these potential 
contributions arise solely from the valence electrons. 

The input potential for most TB calculations is assumed to be the sum of neutral 
atomic potentials (Va + Vb) as in Fig. 1. However, as can be seen from Eq. (38) 
the output potential can differ from this input potential due to the flow of charge. In 
particular, the third contribution in Eq. (38) results from the shift in the potential due 
to the flow of charge q from one atom to another, whereas the fourth contribution 
arises from the creation of the covalent bond with a bond order — q 2 . This latter 
term is given by 

V£ov(r) = J p cov (r')/\r-r'\dr'. (42) 

Thus, for homovalent dimers such as H2 with A = the output potential simplifies 
to 

[V AB (r)]A=o = V A (r) + V B (r) + Vf ov (r) . (43) 

A self-consistent TB calculation requires the output potential to be identical to the 
input potential. This will be discussed in the next section once we have determined 
explicit expressions for the on-site and intersite Hamiltonian matrix elements, and 
hence the appropriate TB parameters. 

The on-site Hamiltonian matrix elements Haa and Hbb follow from Eqs. (7), (9) and 
(38), namely 

Haa = E° A + a A + (J AA - J AB )q + J A cV^ ~ Q 2 , (44) 

where E A is the atomic energy level of the isolated free atom. The second con- 
tribution is the shift in the on-site energy on atom A due to the crystal field from 
neighbouring atom B, that is 

a a = j pAV B dr . (45) 

The third contribution is the additional shift in the on-site energy on atom A due 
to the flow of charge q, which gives rise to an upward shift J AA q due to the in- 
creased number +q of electrons on site A and a downward shift Jab{—i) due to the 
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Coulomb attraction resulting from the negative charge -q on the neighbouring site 
B. 

The on-site Coulomb integral (traditionally denoted by U in the many-body Hubbard 
Hamiltonian) is given by 

Jaa = J J p A (r)pA(r')/\r-r'\drdr' (46) 
and the intersite Coulomb integral by 

J AB (R) = j f P A(r)p B (r' - R)/\r - r'\drdr' . (47) 

The fourth contribution is the Coulomb shift in the on-site energy on atom A due to 
the formation of the covalent bond, where the atom-covalent bond Coulomb integral is 
given by 

Jac = J J PA(r)Pcov(r')/\r-r'\drdr' . (48) 
Similarly the on-site energy on atom B is given by 

H B B = E° B +a B - (J BB - J AB )q + J Bc y/l - q 2 . (49) 

The inter-site Hamiltonian matrix element H AB follows from Eqs. (7), (10) and (38), 
namely 

Hab = Psk + \(E% + E A )S - q j <j> A {V? - V p A )<\> B dr 

- Q f MV& ~ Vfifodr + ^Jl-q 1 j ^V^dr (50) 

where the first term is the well-known Slater-Koster two-centre bond (or hopping) in- 
tegral for orthogonal orbitals, namely 

Psk{R) = J Mr)l(VA + V B )/2]0 B (r-R)dr. (51) 

Because the valence s-orbitals 4>a and 4>b are angularly independent, this bond in- 
tegral is a function only of R, the magnitude of the internuclear distance R, and 
not the direction R. Psk is clearly negative for atomic s orbitals interacting via the 
negative overlap of the atomic potentials (see Fig. 1). 

Finally, using Eq. (31) to express the quantum-mechanical interference factor 4>a{i~) x 
4> B (r — R) in terms of the overlap integral S and the covalent-bond density p cov , the 
intersite Hamiltonian matrix element can be written as 

Hab = /35^+^(^+^)S-^[(J Bc -JA C )+(JBB-JAA)5]g+^[J cc +(J Bc +JAc)^]^l - q 2 , 

(52) 
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where the covalent-bond-covalent-bond Coulomb integral is given by 

J cc = J J Pcov(r)pcov(r')/\r -r'\drdr' . (53) 

The key TB parameters E, A and f3 can now be found by substituting Eqs. (44) and 
(52) into Eqs. (12), (14) and (17) respectively. 

We find the average on-site energy of the dimer is given by 



E = -(H B b + Haa) = E + -(a B + a A ) - ^(Jbb - Jaa)q+ ^{Jbc + Ja c )V^ ~ Q 2 , 

(54) 

whereto = \{E° B + E\). 

On the other hand, the atomic energy-level mismatch in the dimer can be written 

A = H BB - H AA = A + (a B -a A )- 2Jq + {J Bc - J Ac) V^-Q 2 , (55) 
where Ao = (E B — E A ) and the effective Coulomb integral 

J = \( j bb + Jaa) - Jab ■ (56) 
Finally the bond integral can be expressed in the form 

P = H AB -ES= [A, - \{J Bc ~ J Ac )q + \jce\/l-q 2 ] , (57) 

where 

0o = Psk~ \{ctA + a B )S. (58) 

The first term in Eq. (57) is, therefore, the usual Slater-Koster bond integral that has 
been modified by the non-orthogonality of the atomic orbitals. The second term is 
the contribution to the bond integral due to the creation of the point charges ±q , 
whereas the third term arises from the formation of the bond-charge density with 
bond order y/l — q 2 . 



2.4 Self-consistency 

We have seen from the discussion following Eq. (38) that the output potential differs 
from the input potential due to the redistribution or flow of charge reflecting the 
formation of both the covalent bond with bond order — q 2 and the ionic bond 
with point charges ±q. This in turn influences the values of the TB parameters. In 
particular, the atomic energy-level mismatch A depends explicitly on the charge q 
through Eq. (55), 



A(q) = A + (a B - a A ) - 2Jq + ( J Bc - J Ac )^l-q 2 . 



(59) 
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As expected, the right-hand side of this equation involves various Coulomb inte- 
grals that drive the shift in the on-site energy levels. On the other hand, we have 
already found by solving the TB secular equation that A depends explicitly on q 
through the simple relationship of Eq. (34), namely 

A(g) = 2|/?|A = 2|/%/yTV. (60) 

Thus, substituting Eq. (57) for (3 into Eq. (60) and equating the right-hand sides of 
Eqs. (59) and (60), we find the following self-consistency equation for q, namely 

= A + (qg - Oa) + ( Jbc ~ Jac)(1 ~ 2g 2 )/y / l - q 2 

2J-J cc + 2|A)|/v'l-g 2 

Within traditional TB methods the contributions arising from the differences in the 
crystal-fields and atom-bond Coulomb integrals are neglected in the numerator, and 
the small bond-bond Coulomb integral is neglected in the denominator. The simpli- 
fied self-consistency equation then takes the form 



An Ar 



(62) 



2J + 2|/3 |/yr^2 1+21^1/0^2' 

where A and j3 Q have been normalized by 2 J i.e. 

A = A /2J,(3 = (3 /2J. (63) 

This corresponds to the quartic equation 

q 4 = -2A q 3 + [4/3 2 + A 2 - l]q 2 + 2A q - A 2 = . (64) 



Fig. 4 shows the contours of the self-consistent charge q (and corresponding bond 
order 0) as a function of the two variables 2/3n and An , where we have chosen the 
co-ordinates (2/3 , An) in keeping with Eqs. (23)and (62). It follows from Eq. (63) 
that 

Ao = (g/Vl-9 2 )2A) + 9, (65) 

so that the curves in Fig. 4 for different q vary linearly, as observed. The solid curve 
corresponds to q = 6 = l/\/2 when the charge and bond-order take equal values. 
Since the maximum charge and bond-order for an s-valent dimer is 1 we can regard 
q as the degree of ionicity of the bond and 6 as the degree of covalency. Thus, the 
ionicity and covalency increase in the direction of the arrows shown. As expected, 
for homovalent dimers with the normalized atomic energy-level mismatch An = 
we have fully covalent dimers with q = 0, 6 = 1. For the normalized bond integral 
/3n = corresponding to isolated free atoms we have from Eq. (62) that 

for An > 1 
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1 2 

W = Po/J 

Figure 4: Contours of self-consistent charge q (and corresponding bond order = 
a/1 — q 2 ) as a function of 2f3 and A Q . Solid black curve corresponds to q = 6 = l/y/2 
when degree of ionicity equals degree for covalency. Coordinate point for LiH, marked 
by *, corresponds to (1.2, 0.88) in q = 0.4. 

Therefore, the different constant (/-curves in Fig. 4 intercept the vertical axis at 
Aq = q for q < 1. We will return to a discussion of the free atom case and the 
heterovalent dimer LiH in the next section once we have found explicit values for 
the TB parameters Aq , /3o and J. 



3 Binding energy 
3.1 Total energy 

The total energy of the AB s-valent dimer can be written in the form 

U = U band -U dc + l/R, (67) 

where the first contribution is the band energy and the second term is the double- 
counting energy. The third contribution is the core-core repulsive energy (for the 
case of LiH this only holds for internuclear separations R > Rf[ e (c.f. Eq. (40)). We 
have neglected the usual exchange-correlation contributions, assuming that they 
have been absorbed implicitly within the fitted on-site and inter-site TB parameters. 



The band energy is the sum over the occupied eigenvalues, namely 
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Uband — ^2 
n occ 

where the prefactor 2 accounts for the spin-degeneracy in non-magnetic systems 
(c.f. Eqs. (1.175) and (7.13) of [3]). Note that we have retained the nomenclature 
'band' for this sum over occupied eigenvalues even though molecules have a discrete 
spectrum rather than a band of states as in bulk materials. 

For the AB dimer 

U band = 2E+ = 2[\0\S + 1 -{H AA + H BB ) - |/?|/Vl-<7 2 ] , (69) 
where from Eq (33) we have substituted 

^/i + A 2 = l/Vi-g 2 , (70) 

into the eigenvalue in Eq. (19). The first term represents the upward shift in the 
eigenvalue due to the overlap repulsion. The second term locates the absolute energy 
of the average on-site energies with respect to the vacuum level, and the third term 
represents the downward shift in the eigenvalue due to the formation of the covalent 
bond. 

A more transparent understanding of these last two contributions can be obtained 
by rewriting them as 



{HAA + HBB) _ mV ^ . { X^ + -W^7 I ' (71) 

But from Eqs. (14) and (60) 

(H BB - H AA ) = A = 2|/3|A = 2|/%A/l-<7 2 ] • (72) 
Therefore, the band energy can be re-expressed as 

U band = 2\P\S + (1 + q)H AA + (1 - q)H BB - 2\p\y/\-q*\ , (73) 

where from Eq. (33) y/l - q 2 is the bond order G . Thus, the band energy is given by 
the sum of the overlap repulsion, the energy of (1 + q) and (1 - (?) electrons sitting 
on the A and B sites respectively, and the bond energy (given by the product of the 
bond order and the bond integral). 

This physically transparent result for the dimer is a specific example of the well- 
known general result that 

U band = J2 2E n = Kf>H, (74) 

nocc 
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where p and H are the density and Hamiltonian operators respectively (see, for 
example, [3]). Expanding the eigenfunctions in terms of local orbitals as in Eq. (3) 
leads to 

nocc i : j 

assuming real eigenvectors. Thus, for our AB dimer 

U band = 2c\H AA + 2c 2 B H BB + 4c A c B H AB . (76) 
Taking the eigenvectors from Eqs. (21) and (22) we can write 



Uband = f [(1 + q) ~ S ^T^]H AA + [(1 -q)- S^1^]H BB 1 
| +2[ v / l - q 2 - S]H AB J 

where we have used Eqs. (32) and (33). Finally, substituting H AB = (3 + \{H AA + 
Hbb)S from Eq. (17) into the above equation, we recover 

U band = [(1 + q)H AA + (1 - q)H BB - 2^1 - q 2 + 2\(3\S] (78) 

to first order in S, where the repulsive overlap contribution arises from the non- 
orthogonality term in the original bond order 0^ of Eq. (29). Thus, Eq. (78) ob- 
tained from the trace of pH is identical to Eq. (73) obtained from the eigenvalues 
directly, as it must. 

The band energy, Eq. (73), can be expressed explicitly in terms of the TB parameters 
by substituting in both Eq. (54) for the average on-site energy E and also Eq. (55) 
for the atomic energy-level mismatch A. It follows that 



= < 



U band = 2\f3\S + {H BB + H AA ) - q(H BB - H AA ) - 2|/3|^I - q 2 (79) 
2\P\S + (E% + E° A ) + (a B + a A ) - {J BB - J AA )q \ 
+ {Jbc + J Ac )^fl -q 2 - q[A + (a B ~ a A )\ 
-2[\^\ + \{a A + a B )Syi-q 2 
+J cc (l-q 2 ) + 2Jq 2 

This band energy has double-counted the Coulomb energy between the electrons 
so that the total energy in Eq. (67) includes the double-counting correction term, 
namely 



Udc = \j j p(r)p{r')/\r-r'\drdr' . 



(80) 



Substituting in the valence charge density from Eq. (30), this can be written in terms 
of the various Coulomb integrals as 



Udc 



\J AA {\ + q) 2 + \{l - q) 2 J BB + (1 - q 2 ) J AB 

+v / i-<? 2 [(i + q)JAc + (i - q)J Bc ] + |(i - q 2 )Jc 



(81) 
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The terms in the above double-counting expression can be regrouped to give 

{[\{Jbb + Jaa) + Jab] + \J C c 1 
- {Jbb - Jaa)q - {J Be + JaJV 1 - Q 2 } >• (82) 
+ ( Jbc ~ JaM^I^ + ( J ~ \Jcc)q 2 J 

The first contribution is the electronic Coulomb energy of neutral atoms A and B as- 
suming the atomic charge density remains frozen as the atoms are brought together 
from infinity to form the AB dimer. The second contribution is the Coulomb self- 
energy of a neutral covalent bond that is fully saturated with a bond order of unity. 
The next two terms in the square brackets reflect the first order shifts in the average 
on-site energy in Eq. (54) due to the charge q and bond order = \Jl — q 2 cor- 
recting the dimer potential in Eq. (38) beyond that of overlapping neutral frozen 
atoms. The following contribution is a cross-term that arises from the output charge 
q interacting with the covalent-bond potential in Eq. (38), and vice versa. The last 
contribution in Eq. (82) is the second-order q 2 double-counting Coulomb terms. 

The total energy of the AB dimer, Eq. (67), then takes the physically transparent 
form 



U(R, q) = [(U° B + U° A ) + U over (R) + U es (R) + U cov (R, q) + U ionic (R, q)] , (83) 
where (Ug + U^) is the energy of isolated neutral free atoms A and B, namely 



(84) 



U°a = E\ - \J AA 

U B = E B ~ 2 J BB 

The next four contributions on the right-hand side of Eq. (83) are the overlap, elec- 
trostatic, covalent and ionic terms respectively. They are given to first order in the 
overlap S by 

U mer (R) = 2\Po\S, (85) 
U es (R) = (1/R + a B + a A - J AB ) + - J cc , (86) 
U cov (R,q) = -2\(3 \^/l-q 2 , (87) 

and 

U lomc (R, q) = {- [A ° + ^ B ~ ( /t 2 X (88) 

I +iy 2 'Jccjq ) 

All four contributions are functions of the internuclear separation R since all the TB 
parameters except Ao on the right-hand sides are distance dependent. The covalent 
and ionic contributions are also explicitly dependent on the charge q. 



We see that the non-orthogonality of the atomic orbitals enters only the repulsive 
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overlap and attractive covalent-bond contributions explicitly. The former represents 
the well-known quantum mechanical overlap repulsion, whereas /3q in the latter is 
the usual Slater-Koster bond-integral modified by the shift \ (olb + cya)S through Eq. 
(58). The origin of both these overlap terms can be understood by considering the 
homovalent s-valent dimer with q = 0. In this case, the non-orthogonal TB secular 
Eq. (8) has matrix elements 

Haa = H BB = E° + a (89) 

from Eq. (44) and 

Hab = Psk + E°S (90) 

from Eq. (52), where we have made the common TB assumption that Vab = Va + Vb 
by neglecting the small covalent-bond contribution to the dimer potential. 

These matrix elements lead to the eigenvalues 

E ± = E° + (a ± p SK )/(l ± S) . (91) 
These can be written in the more familiar form as 

E ± = E° + a* ± /3 (92) 
by defining an effective crystal field term 

a* = \{E + + E~) -E° = (a- P SK S)/{1 - S 2 ) (93) 
and bond integral 

P* = \{E + - E-) = {Psk - aS)/(l - S 2 ) . (94) 

It follows from Eq. (58) that a* = (a — floS) and /3* = (3q to first order in S. Thus, we 
see that within the total energy expression, Eq. (83), the shift in the eigenvalues a* 
enters the overlap repulsion through 21/3015" in Eq. (85) and the electrostatic contri- 
bution through («b + (xa) in Eq. (86) (for the 2 valence electrons per dimer). On the 
other hand, the bond integral (3q enters the covalent bond contribution as 2|/3o|0- 

The electrostatic contribution, Eq. (86), comprises the Coulomb energy of overlap- 
ping frozen neutral atoms (the bracketed term) and the Coulomb self-energy of the 
neutral covalent bond. This is complemented by the ionic contribution, Eq. (85), that 
accounts for the flow of charge q that is initially driven by the atomic energy-level 
mismatch Ao- Whereas this leads to the first-order contribution — Aog in Equation 
(88), the other first-order contributions —{Jbb — Jaa)q and —(Jbc — J Ac) \A — Q 2 
in the band energy, Eq. (79), are cancelled by corresponding terms in the double- 
counting energy, Eq. (82). These had arisen in the band contribution from the 
shift in the average on-site energy in Eq. (54) due to the charge q and bond or- 
der 6 = y/l — q 2 correcting the dimer potential beyond that of overlapping neutral 
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frozen atoms. This important result is an example of the Harris-Foulkes DFT func- 
tional that is accurate to first order even though its input potential is that due to 
overlapping frozen atoms (c.f. section 3.4 of [33]). 

For a given internuclear separation the total energy will be stationary for the value 
of q that satisfies 

d MM = 0. (95) 
oq 

It is trivial to show that this stationary condition for the total energy leads to the 
same self-consistency equation for q that was obtained by solving the TB secu- 
lar equation self-consistently namely Eq. (61). This confirms that we have in- 
deed solved the AB s-valent dimer problem correctly within our assumptions of 
a minimal basis and neglect of explicit exchange-correlation terms. Just as this self- 
consistency equation, Eq. (61), was simplified to Eq. (62) by making well-justified 
approximations, so traditional TB approximates Eqs. (85) - (88) by neglecting the 
small contributions arising from the differences in the crystal fields, (erf - a^) and 
the Coulomb integrals involving the covalent-bond density p cov . That is, the total 
energy is given by Eq. (83) with the electrostatic, covalent and ionic contributions 
approximated by 

U es {R) = l/R + a B + a A - J AB , (96) 
U cov (R,q) = -2\p \Vl-q 2 , (97) 

and 

U iomc (R,q) = -A q + Jq 2 . (98) 



Finally, the force is given by 



F _ dU[R,g(R)} (99) 
dR 

dU dU dq _ _dU 
~ ~dR ~ ~dqd~R ~ 



because the second term vanishes by Eq. (95). Thus, 

-[U' over (R) + U' es (R)} 
-2\p' {R)\y/T=tf - j' AB (R)q 



2 



(100) 



where g = dg/dR and we have used the approximate expressions (96) - (98). This is 
consistent with the Hellmann-Feynman theorem in that the force does not depend 
on the derivative of the density matrix elements (sections 3.1 and 7.6 of [3]). From 
Eqs. (76) - (78) this is equivalent to the force being independent of either the first- 
order change in the charge q or the bond order = y 7 ! — q 2 . 
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3.2 Free-atom limit 

In the limit as the internuclear separation tends to infinity the total energy in Eq (83) 
tends to the free-atom limit 

U(R -> oo, q) = {U A + U° B ) - A q + ^(J BB + J AA )q 2 . (101) 

This corresponds to the energy of free atoms A and B with net charges Qa = —q and 
Qb = q respectively where 

Ua{Q) = U a + \E a \Q + \.JaaQ 2 \ nm . 
U B {Q) = U° b + \E%\Q + \J BB Q 2 J • KWA) 

This second-order expansion of the free-atom energy in terms of the net charge Q 
allows the TB linear prefactor | E A , B s | and quadratic prefactor \ Jaa( b b) to be identi- 
fied with the experimental ionization potential (IP) and electron affinity (EA) of the 
respective atoms. The latter are defined by 

IP = U(Q = 1) - U{Q = 0) (103) 

and 

EA = U(Q = 0) - U(Q = -1) , (104) 
where U(Q) is the energy of the charged atom. It follows from Eqs. (102) - (104) that 

\E°\ = ^IP + EA) = x° (105) 

and 

l -J*= l -(IP-EA) = r } (106) 

with J° = Jaa(bb)- X° is the Mulliken electronegativity which he defined as the 
average of the ionization potential and electron affinity of the free atoms (page 107 
[43]). rj is called the chemical hardness of the free atom and is defined by the quadratic 
prefactor in the free-atom energy expansion, Eq. (102) (page 107 [4]). Thus, within 
this simple s valent TB model the magnitude of the energy level of the neutral free 
atom plays the role of electronegativity, whereas one-half of the on-site Coulomb 
integral (or Hubbard U) plays the role of hardness. 

The change in total energy with the flow of charge can be written from Eq. (101) in 
the free-atom limit as 

AU(R^ oo,q) = -(Ax )q + 2r]q 2 (107) 

where Ax° = x% ~ Xa anc ^ V = \ (vb + Va) ■ The values of x° and r\ for H and Li are 
given in table 1. 
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IP (eV) 


EA (eV) 


x ° = E° = \{IP + EA) 


t, = ±J°=±(IP-EA) 


H 


13.6 


0.8 


7.2 


6.4 


Li 


5.4 


0.6 


3.0 


2.4 



Table 1: Experimental values of ionization potential, electron affinity, and correspond- 
ing values of electronegativity (magnitude of free-atom energy level) and chemical 
hardness (one-half on-site Coulomb integral) for hydrogen and lithium. 



Fig. 5 shows the resultant energy curve AU (R — > oo, q) for LiH with Ax° = 4.2 eV 
and fj = 4.4 eV at infinite internuclear separation. We see that charge is driven from 
the Li atom to the H atom by the attractive linear term proportional to their elec- 
tronegativity difference Ax° , but that this is countered by the repulsive quadratic 
term proportional to their average chemical hardness fj. Provided Ax°/4fy < 1 the 
total energy Eq. (107) is stationary for 

A X ° A 

Qoc = -Tzr = -j —j — , (108) 

4ry Jaa + Jbb 

when the binding energy with respect to the neutral free atoms takes the Pauling 
form for the heat of formation [5], namely 

AH = -l(A x ) 2 . (109) 

These last two expressions give q^ = 0.24 and AH = -0.50 eV for the case of LiH 
in the free-atom limit. In reality when the Li and H atoms are separated to infinity 
the isolated atoms cannot contain fractional charges since there will be no hopping 
or tunnelling of the electrons between the sites. The only physical configurations 
would be Li°H°, Li + H~ and Li~H + . But from Fig. 5 we see that the configura- 
tion Li + H~ costs nearly 5eV in energy compared to the charge neutral state, and 
Li~H + would still be a further 8.4eV higher in energy. Thus, energetically nature 
favours the separation of the LiH dimer into neutral atoms. The prediction of frac- 
tional charges in the free-atom limit is a well-recognized problem with TB and DFT 
because they have not imposed any constraint on integer numbers of electrons. We 
should, therefore, be aware that these one-electron-type methods can lead to an 
over-estimation of charge flow at large internuclear separations (Eq. (46) of [4]). 



3.3 Binding energy 

The binding energy of the AB dimer is obtained by subtracting off the energy of the 
neutral free atoms A and B, Eq. (84), from the total energy, Eq. (83). That is, 

U be (R, q) = U rep (R) + U C0V (R, q) + U ionic (R, q) , (110) 
where the repulsive contribution is given by 
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Figure 5: Changes in total energy of LiH in limit of infinite internuclear separation as 
function of net charge q on Li atom. A% = {E^ — E° H | is electronegativity difference and 
77 = \{Juia + Jhh) is average chemical hardness. Open square marks position where 
energy is stationary at = 0.24. 



U rep (R) = U over (R) + U es {R) . (Ill) 
We will make the customary TB approximation as in Eqs. (96) - (98) to write 

U cov (R,q) = -2\(3 (R)\ v / l-q 2 , (112) 

and 

U ionic (R, q) = -A q + J{R)q 2 . (113) 

As we have seen in the previous section this binding-energy expression tends to 
the wrong limit, Eq. (108), as the atoms are pulled infinitely apart. However, this 
need not be of a major concern to us as TB and DFT are known to give an accurate 
representation of the energetics of most materials in the vicinity of their equilibrium 
separations or atomic volumes, where the effective one-electron-type description 
works well. 

Let us consider first the homovalent hydrogen molecule since the ionic contribution 
vanishes as Ao = 0. The upper panel of Fig. 6 shows the bonding and anti-bonding 
energy levels that are predicted by non-spin polarized DFT. This neglects the occur- 
rence of local magnetic moments that are stable for internuclear separations greater 
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Figure 6: TB fit to DFT energy levels and binding energy curve as function of internu- 
clear separation R for H 2 with equilibrium distance marked by dashed vertical line [6]. 
Upper panel: DFT bonding (E + ) and anti-bonding (E~) energy levels. Middle panel: 
TB parameters E°, a* and (3q from fitting DFT data in upper panel. Lower panel: DFT 
binding energy (solid black curve), covalent-bond energy — 2|/3 | (solid blue curve), re- 
pulsive energy (solid red curve) with its overlap (dotted red curve) and electrostatic 
(dashed red curve) components. 



than about 2R$, where Ro = 1.4 au is the equilibrium distance (see, for example, Fig. 
3.6 of [1]). 

These two DFT eigenvalues can be used to fit the three TB parameters E°, a* (R) and 
Po{R) because both a*(R) and /3o(i?) tend to zero at large internuclear separations, 
as shown in the middle panel of Fig. 6. We see that the repulsive overlap contribu- 
tion \f3o\S in Eq. (93) dominates over the attractive crystal field term a in Eq. (92) as 
a* > over the range of distances plotted. However, we observe that for R < R^a* 
starts to saturate as the third order attractive contribution a* S 2 in Eq. (93) becomes 
important for S > ^. At the equilibrium seperation a* = 2.0eV and /3o = 5.8eV. 

Having used the DFT energy levels to fit the TB parameters E°, a* and /3q, we can 
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now use the DFT binding energy curve in the lower panel of Fig. 6 to extract the 
repulsive energy, from which we may deduce the corresponding overlap and elec- 
trostatic contributions. The repulsive energy in Eq. (Ill) can be expressed as 

U rep (R) = U be (R)-U cov (R), (114) 

where U cov (R) = 2\(3 (R)\ from Eq. (112) as the hydrogen molecule has a fully 
saturated bond with a bond order of unity. Subtracting the resultant covalent-bond 
contribution in Fig. 6 from the binding energy we recover the repulsive curve that 
is plotted. 

This repulsive contribution may be separated into its overlap and electrostatic con- 
stituents as follows. Firstly, the overlap term can be deduced by assuming that the 
overlap integral is approximately proportional to the bond integral, so that the over- 
lap repulsion varies as 

U over (R) = A[f3 (R)] x , (115) 

The dotted curve in Fig. 6 demonstrates an extremely good fit to the DFT-derived 
repulsion for R > Rq with A = 0.185 and A = 1.85 (giving the overlap energy in Eq. 
(115) in units of eV/dimer provided (3q is in eV). This value of 1.85 is remarkably 
close to the original Wolfsberg-Helmholz approximation of A =2 [1]. Secondly, the 
electrostatic contribution can be obtained by subtracting the overlap repulsion from 
the DFT-deduced repulsive curve in Fig. 6. The resultant dashed curve kicks-in for 
R < Ro and has a shape that agrees with that obtained analytically by evaluating 
the electrostatic interaction between neutral atoms comprising Is atomic charges 
resulting from variational hydrogonic wave functions (C 3 ) 5 exp (-(r) (c.f. Eq. (32) 
and Fig. 6 of [7]). Traditional TB methods write the repulsive energy as a sum over 
pair potentials, namely 

U rep (R) = ±J2<S> rep (R). (116) 

That is, for the dimer U rep (R) = § rep (R) This pairwise assumption is justified to 
lowest order in bulk materials: the electrostatic interaction is pairwise within the 
usual TB approximation that the bulk potential is the sum over atomic potentials 
and the overlap contribution, resulting from the non-orthogonality shifts in the on- 
site energies, is also pairwise as Yli j HijSji to first order in S. 

We see, therefore, that for homovalent dimers the binding energy can be decomposed 
in terms of the physically-based contributions of overlap repulsion, atom-atom elec- 
trostatic interaction, and covalent bond energy. The firm physical basis of this TB 
model is confirmed by the fact that the bond integrals extracted from the DFT eigen- 
values leads to a covalent-bond contribution to the energy that reproduces the DFT 
binding energy curves almost exactly for R > 2R , when at smaller distances the 
repulsive overlap energy becomes noticeable. Further, the electrostatic atom-atom 
interaction is found to provide the rapid increase in repulsion at short distances. 
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This soundness of the TB model underpins its application to many material systems 
where it is assumed that the two-centre bond integrals Po(R) and the pairwise re- 
pulsion § rep (R) are transferable from one crystal structure or atomic configuration 
to another. 

Finally we consider the heterovalent dimer LiH where the ionic contribution also 
enters the TB expression for the binding energy Eq. (110). We have seen from Eq. 
(62) that the self-consistent charge q is a function of the normalized bond integral 
A) = /3o/2J and atomic energy-level mismatch A = A /2J , where J = \{Jbb + 
Jaa) - Jab(R)- From table 1 the values of A = E% - E Q A = 4.2eV and \{Jhh + 
JuLi) = 8.8eV. The value of the intersite Coulomb integral Jab(R) can be found 
from the interpolation formula [8] 

Jab(R) = 1 (117) 

[R 3 + (1/VJaaJbb) 3 P 

using atomic units with R in au and the Coulomb integrals in Hartrees. For large 
distances when the A and B atoms do not overlap this falls off Coulombically as 
1/R. However, for short distances it saturates as the atomic orbitals overlap. For R 
= it takes the value V JaaJbb which is exact for the homovalent case when the A 
and B atoms are chemically identical. From table 1 for LiH V 'JaaJbb = 7.8eV= 0.29 
Hartree. Thus, 

27 2 

Juh(R) = : r eV . (118) 

[#3 + 41.9] a 

This is plotted in Fig. 7 where we see that the saturation from the unscreened 
Coulomb inverse power law becomes very important for R < 5 au. 

Table 2 gives the DFT values of the equilibrium internuclear separation Rq and bond 
integral /?o for H2 and Li2. The equilibrium separation for LiH was taken from 
Zen's law as the arithmetic mean of the equilibrium separations for Li2 and H2. 
The equilibrium bond integral for LiH, on the other hand, was assumed to be the 
geometric mean of the equilibrium bond integrals for Li2 and H2 (as for the LiH on- 
site Coulomb integral in Eq. (117)). We see from Fig. 7 that the equilibrium value 
of J LiH is then 6.4eV. Hence, J in Eq. (62) takes the value (8.8 - 6.4) = 2.4eV, so that 
the presence of the attractive nearest neighbour Coulomb interaction dramatically 
reduces the effective on-site Coulomb integral, thereby enhancing the propensity for 
charge flow. Thus, LiH at equilibrium takes the co-ordinate point (2/3q, Aq) = (1.20, 
0.88) on the charge contour map in Fig. 4. This corresponds to a value of q = 0.4, so 
that, as expected, the ionic and covalent bonds in LiH play an almost equal role. 



4 Conclusions 



These notes have shown how it is possible to derive rigorously a TB expression 
for the binding energy of an AB s-valent dimer by starting from the effective one- 
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5 10 

R(au) 

Figure 7: Comparison of LiH interatomic Coulomb integral Juh(R) with unscreened 
Coulomb inverse power law for a wide range of internuclear seperations R. Vertical 
dashed line gives equilibrium internuclear separation. 





H 2 


Li 2 


LiH 


-Ro(au) 


1.4 


5.1 


3.3 




5.8 


1.5 


2.9 



Table 2: Equilibrium bond lengths and bond integrals for H 2 , Li 2 and LiH. 



electron-type equations of DFT. The TB binding energy has four well-defined contri- 
butions that are physically and chemically motivated, namely the overlap repulsion, 
the atom-atom electrostatic interaction, and the attractive covalent and ionic bond 
energies respectively. This TB binding energy in its turn can then be used to justify 
well-founded ionic or covalent interatomic potentials. Hopefully these notes will 
help you the reader to hone your own physical intuition about the nature of chemi- 
cal bonding in materials, and hence to choose your own path through the myriad of 
different TB approaches out there in the literature. 
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